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Abstract 

Theoretical and numerical studies are performed for the nonlinear structures (explosive, solitons 
and shock) in quantum electron-positron-ion magnetoplasmas. For this purpose, the reductive per- 
turbation method is employed to the quantum hydrodynamical equations and the Poisson equation, 
obtaining extended quantum Zakharov-Kuznetsov equation. The latter has been solved using the 
generalized expansion method to obtain a set of analytical solutions, which reflect the possibility 
of the propagation of various nonlinear structures. The relevance of the present investigation to 
the white dwarfs is highlighted. 
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I. INTRODUCTION 

Numerous investigations relating to wave phenomena, have been studied 

in dense quantum plasmas, are of fundamental importance for understanding collective in- 
teractions in superdense astrophysical environments {(J, in high intense laser-solid density 
experiments 7], in ultracold plasmas ]8], in microplasmas and in micro-electronic de- 
vices [h]]. New characteristics of quantum plasma arise due to the pressure law describing 
the fermionic behavior of the charged carriers, quantum forces associated with the electron 
tunneling, as well as the Bohr magnetization involving the electron 1/2 spin. The quantum 
Bohm potential produces modifications in the dispersions of collective modes at quantum 
scales. The latter are strongly effected by the plasma number densities and Fermi tempera- 
tures. It is well-known that quantum mechanical effects become relevant when the thermal 
de Broglie wavelength of the charged particles is equal or larger than the average interpar- 
ticle distance. In particular, quantum behavior of the electrons reaches much easily due to 
less mass compared to ions. 
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15| have been 



In recent years, many theoretical and numerical analysis jll . 
carried out to investigating the new features of plasmas with quantum corrections by using 
both the Schrodinger-Poisson and the Wigner-Poisson systems. In this context, Manfredi 11] 
reported different approaches to model the collisionless electrostatic dense quantum plasmas. 
Haas et al. [12| investigated the linear and nonlinear properties of the quantum ion-acoustic 
(QIA) waves in dense quantum plasmas by employing the quantum hydrodynamical (QHD) 
equations for inertialess electrons and mobile ions. They examined that the quantum Bohm 
potential modifies the linear wave dispersion and affects strongly the QIA solitary waves. 
Shukla and Eliasson [3] presented the numerical study of the dark solitons and vortices in 
quantum electron plasmas. Moslem et al. investigated the quantum dust-acoustic double 
layers in a multi-species quantum dusty plasma. It was found that both compressive and 
rarefactive double layers can only exist for positively charged dust particles. Later, Ali et al. 
[if] ] studied the QIA waves in a three-component plasma, comprised of electrons, positrons, 
and ions. They employed the reductive perturbation method and pseudo-potential approach 
for the small and arbitrary amplitude nonlinear QIA waves, respectively. It was shown that 
the amplitude and width are significantly altered due to the quantum statistics and quantum 
tunneling effects. Misra et al. 17j considered the nonlinear propagation of electron-acoustic 
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waves in a nonplanar quantum plasma, consisting of two groups of electrons: the inertial 
cold electrons and inertialess hot electrons as well as the stationary ions. They obtained the 
bright and dark solitons depending strongly upon the presence of cold electrons. 

The laboratory and dense astrophysical quantum plasmas can be confined by an external 
magnetic field. Therefore, the effect of the magnetic field has to taken into account, espe- 
cially for astrophysical observations (such as white dwarfs, neutron stars, magnetars, etc.) 
where the high magnetic field plays an important role in the formation and stability of the 
existing waves. Several authors have considered the effect of magnetic field in different quan- 

n 

turn plasma models. For example, Haas [18j introduced a three-dimensional QHD model 
for dense magnetoplasmas and established the conditions for an equilibrium in the ideal 
quantum magnetohydrodynamics (QMHD). Ali et al. [3] employed the QMHD equations 
presenting a fully nonlinear theory for ion-sound waves in a dense Fermi magnetoplasma. 
' t was revealed that only subsonic ion-sound solitary waves may exist. Shukla and Stenflo 
2(3] derived the dispersive shear Alfven waves in a quantum magnetoplasma, incorporat- 
ing the strong electron and positron density fluctuations. The shear Alfven modes acquire 
additional dispersion due to quantum corrections. Later, Ali et al. [21[ have been carried 
out for the low-frequency electrostatic drift-like waves in a nonuniform collisional quantum 
magnetoplasma. It was shown that the modes become unstable and can cause cross-field 
anomalous ion- diffusion. 

Three decades ago, Zakharov and Kuznetsov 22| derived an equation for nonlinear ion- 
acoustic waves in a magnetized plasma containing cold ions and hot isothermal electrons. 
The Zakharov-Kuznetsov (ZK) equation has also been derived for different physical systems 



and scenarios 
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24j . Nonlinear wave solution for Z 



^ equation can produce an instability 



in a three-dimensional system as discussed in Refs. [25l |26|. Moslem et al. 27| extended 
the work for a three-dimensional nonlinear ion-acoustic waves in a quantum magnetoplasma, 
highlighting the bending instability of the solitary wave solution of the quantum ZK equa- 
tion. Recently, Masood and Mushtaq [28] studied obliquely propagating electron-acoustic 
waves in a two-electron population quantum magnetoplasma and examining the effects of 
nonlinearity at quantum scales. 

In the present paper, we shall investigate the possible nonlinear structures (soliton, ex- 
plosive and shock pulses) of the QIA waves in a collisionless electron-positron-ion magneto- 
plasma using the QHD equations. By means of computational investigations, we examine 
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the effect of the positron concentration, the quantum diffraction and the quantum statistical 
effects on the profiles of the nonlinear excitations. The paper is organized as follows: The 
basic equations governing the dynamics of the nonlinear Q1A waves are presented and the 
extended quantum ZK equation describing the system is derived in Sec II. In Sections III 
and IV, we apply the generalized expansion method to solve the extended quantum ZK 
equation. A set of analytical solutions is obtained, and then used to investigate numerically 
the effect of positrons and the quantum parameters on the nonlinear excitations. The results 
are summarized in section V. 



II. BASIC EQUATIONS AND DERIVATION OF THE EXTENDED QUANTUM 
ZK EQUATION 

We consider a dense magnetoplasma whose constituents are the electrons, positrons, and 
singly charged positive ions. The plasma is confined in an external magnetic field H = H z, 
where z is the unit vector along the z— axis and H is the strength of the magnetic field. 
We assume that the quantum plasma satisfies the condition T Fep ^> T Fi , and obeys the 
electron/positron pressure law P e , p = m V Fe /2>n 2 e ^ where Vpe-p = (ZKbTfcp/ 'M) 1 ^ 2 is 
the electron/positron Fermi thermal speed, Kb is the Boltzmann constant, Tp e ,p {Tn) is 
the electron/positron (ion) Fermi temperature, M is the electron and positron mass, n ejP 
is the electron/positron number density, with the equilibrium value n eiP o- The nonlinear 
propagation of the QIA waves is governed by the dimensionless hydrodynamics equations as 

^ + V. (n t uO = 0, (1) 

^ + Ui.Vui = -V0 + Ui x z, (2) 

at 

Q V 2 = n e - n p - rii, (3) 

n e = /i e (l + 20 + tf e 2 ^py, (4) 



and 

l-2a<P + aH 2 e J , (5) 

where n^, Uj, and are the ion number density, the ion fluid velocity, and the electrostatic 
potential, respectively. Since, the ion mass is much larger than the electron/positron mass, 



one can ignore the quantum effects of the ions in Eq. (2). The statistical and diffraction 
effect for the system can be seen through the nondimensional parameters a(= Tpe/Tpp) and 
H e (= eH h/2cy r MiMKBTFe), respectively, where h is the Planck constant divided 
by 2ir, Mi (M) is the ion (electron/positron) mass, and c is the speed of light 
in vacuum. Here, Q(— uj c i/uj p i), where uj a i (= eHo/rriiC) and u P i (= ^J^nehi^jMi) 
are the ion gyrofrequency and the ion plasma frequency, respectively, riio is the 
equilibrium ion density. Equations (4) and (5) reveal that the electrons and positrons 
do not follow the Boltzmann law contrary to the classical plasma. The physical quantities 
appearing in Eqs. (l)-(5) have been appropriately normalized: fi e ^^ p — * ^e^p/^iOj Uj — > 
Uj/C s , t — > tu) ci , V — > Vp s , and — > e<p/2K B T Fe , where p s (= C s /uj ci ) is the ion-sound 
Fermi gyroradius and C s {= ^/2KBTp e /Mi) is the ion-sound Fermi speed. 

Before going to the nonlinear developments, it is necessary to examine the condition for 
neglecting the source term in the continuity equation due to annihilation of plasma species. 
The details are given in the Appendix. 

To investigate the propagation of QIA waves, we expand the dependent variables n e ,i,p, 
Ui, and (j) about their equilibrium values in power of e, 



m = l + en ix + e 2 n i2 + e 3 n i3 + 
n e , P = He, P + en 6iP i + e 2 n e , p2 + e 3 n e , p3 + 



u 



ix, y 



2 3 4 

^ ^ix,yl ^ ^ix,y2 € ^ix,y3 



u iz = eu izl + e 2 u iz2 + e 3 u iz3 + ... 



e0i + e 2 4> 2 + e 3 4> 3 + 



(6) 



where e is a keeping order parameter proportional to the amplitude of the perturbation. 
Following the reductive perturbation method 29], we express the independent variables into 
a moving frame in which the nonlinear wave moves at a phase-speed of A (normalized with 
the ion-sound Fermi speed C s ) as 



X = ex, Y = ey, Z = e(z- Xt) and T = eH. (7) 

The neutrality condition at equilibrium reads /i e = 1 + /i p , where /i e = n e0 /riio and \i v = 
n p o/n i0 . Subistituting (6) and (7) into Eqs. (l)-(5), we obtain the lowest-order in e as 
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1 dfa 
Uiyl = -q^, Uizl = J<Pl, (8) 



n el = // e 0i, n p i = -a/ipfa, 

along with the phase speed rule 

1/2 



It is clear here that the phase speed A of the QIA waves is affected by the quantum statistical 
effect and by the positron concentration /i p . To the next-order in e, we have 



n i2 = 77TT01 + T7T02, U ix2 = A- 



3A 4Y1 A 2 <9X<9Z <9Y" 

M - 2 = + ax ' M " 2 = 2A^ + a (10) 

n e2 = (0 2 - 20 2 ) , n p2 = (a0 2 + 20 2 ) , 



while the Poisson equation gives 



where 



Q0 2 = O, (II) 



[(a 2 -l)^A 4 -A 4 -3] 
W _ 2A* 

Since 0i 7^ 0, therefore Q should be at least of the order of e. Therefore, Q4>\ becomes of the 
order of e 3 ; so it should be included in the next order of the Poisson equation. The next- 
order in e gives a system of equations. Solving this system with the aid of Eqs. (8)-(10), we 
finally obtain the extended quantum ZK equation as 

| +( ^ W) g +c |V +D ^ + |^ 0i (12) 

where we have replaced <pi by <p for simplicity. The nonlinear and dispersion coefficients are 
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given, as 

A 4 + 3 - (<r 2 - 1) /^A 4 



A 
B 



2A 

3((a 3 + 1) /i p A 6 + A 6 - 5) 



4A 3 

C = ^A 3 (40 - (fi P + 1) Hi - a 2 H 2 ^ p ) , 
D = C + ^\ 3 

The extended quantum ZK equation (12) constitutes the final outcome of this model. The 
anticipated balance between dispersion and nonlinearity (which contain the quantum me- 
chanical effects) within the extended quantum ZK equation may give rise to different non- 
linear structures. Some of these solutions will recover in the next section. 

III. EXACT SOLUTIONS OF THE EXTENDED QUANTUM ZK EQUATION 

To obtain the possible analytical solutions of Eq. (12), we assume that 

£ = L X X + L Y Y + L Z Z — #T, (13) 

where Lx, Ly and Lz are the direction cosines and $ is the QIA wave speed to be determined 
later. Using (13) into (12), we obtain 

- #<p' + A <p ip' + B tp 2 ip' + = 0, (14) 
where A = ALz, B n = BL Z and 7 = CL% + DLz (L x + L Y ). According to the generalized 



expansion method 30( the solution of Eq. (14) can represent by 

ip = a + clxco, (15) 

with 

— = k [Co + Ciuo + c 2 uj + C3UJ + c 4 co> ) , (lb) 

where a , 01, c , Ci, c 2 , c 3 and c 4 are arbitrary constants to be determined later and k — ±1. 
Substituting Eq. (15) into Eq. (14) and making use of Eq. (16), we obtain a polynomial 
equation in u. Equating the coefficients of different powers of u, we obtain an overdetermined 
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system of algebraic equations which can be solved with the help of symbolic manipulation 
package Mathematica to give three Jacobi elliptic doubly periodic type solutions as 



A ° _l u I 6 7 c 2 m 2 / / c 2 ~ 



Com 2 (1 — m 2 ) , . 

with c = —2—^ -r , c 2 > 0, c 4 < 0, (17) 

c 4 (2m^ — 1) 



c 2 Cl - m 2 ) 

with c = rA, c 2 > 0, c 4 < 0, (18) 

c 4 (2 — m z ) 



and 



. , / 6 7 c 2 m 2 / / - c 2 A 

Til ^ 

with c = ——2- -2, c 2 < 0, c 4 > 0, (19) 

c 4 (to 2 + 1) 

where m is a modulus of the Jacobian elliptic function and c\ = c 3 = 0. As to — > 1, the 
Jacobi doubly periodic solutions (17) and (18) degenerate to the bell-shapped solitary wave 



^ = "2l| W^bT sech ^O, (20) 
where the arbitrary constant c vanishes. Again, as to — > 1 the solution (19) can degenerate 
to the kink-type wave solution 



^-stW^Hv^ 1 - (21) 

where c = c 2 ,/^. In the solutions (17)-(21), the QIA wave speed -d — \ (— Al/2B + 27c 2 ) 
where c 2 ^ A^/A^Bq. 

Furthermore, the generalized expansion method provides us with further analytical solu- 
tions of the extended quantum ZK equation (12) as 

2 c 2 



c 3 + k a/c 2 - 4c 2 c 4 cosh (2 v /c 2 ~£) 



7? A £? 



with c = ci = 0, c 2 = -, c 3 = - — , c 4 = - — , (22) 

7 07 07 
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and 



2B 



1 + k coth 



-A 



24jB, 



with $ 



42 



and -B < 0. 

(23) 



IV. PARAMETRIC ANALYSIS FOR WHITE DWARFS 

It is clear that the propagation speed of the QIA wave is modified by the effect of the 
quantum statistical effect a and by the presence of positrons fi p . As o and \i p increase, the 
propagation speed of the QIA wave will decrease. The dependence of the nonlinear structures 
amplitude and width on the equilibrium positron number density (fi p ) and quantum effects 
a and H e is more perplex. First, it is important to note that changing [i p leads to a change 
in the phase-speed (A) of the QIA waves [see Eq. (9)] , as well as the electron concentration 
(via the charge-neutrality condition fi e = 1+ fi p ). Since the electron (positron) Fermi 
temperature depends upon the equilibrium electron (positron) number density, it can also be 
affected by ji p through the charge- neutrality condition. As a result, the quantum statistical 
(a) and diffraction (H e ) effects will vary with the positron concentration fi p . 

Based upon the above findings, we shall now investigate the effects of the relevant phys- 
ical quantities, namely the positron concentration n P on the profiles of the QIA nonlinear 
structures. We have used, as a starting point, a typical set of plasma parameter values 



for white dwarfs 



ll| (in the absence of positrons; fi p = 0), namely: n e0 = 4 x 10 28 cm 



4.9 x 10 s K, u ci 



1.88 x 10 ib s- L and u. 



pi 



2.63 x 10 17 s 



17 „-l 



However, once the 



positrons species density is determined, the values of Tp e , A and H e are subsequently com- 
puted, according to the above formulae, which also determine A, B, C and D. In the plots, 
we shall change the positrons concentration, which leads to recalculate all the physical pa- 
rameters again. Obviously, by varying the positron concentration, we simultaneously modify 
all the parameter values used in the plots below. 



A. Solitary and Explosive/Blowup Excitations 

It may be appropriate to point out that the analytical solutions in Sec. Ill have been 
obtained for different arbitrary constants k, Co, ...C4. One of them is the localized solution 
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(22), which is a bell-shapped solitary wave solution. Recall that the arbitrary constant 
k can be either +1 or — 1. For k = — 1, a positive solitary pulse can propagate and for 
k = +1, a negative solitary pulse exist. Note that we have executed the negative solitary 
pulse since it is not physically correct in the the present model. Figure 2 depicts the QIA 
solitary pulse for different values of positron concentration /i p , which now determines Tp e p 
(and the ratio a) through the charge-neutrality condition. It is found that the amplitude of 
the soliton pulse decreases by increasing resulting an increase (decrease) of the electron 
Fermi temperature T Fe (quantum diffraction effect H e ). Physically, the increase of T Fe leads 
to an increase of the electron Fermi energy (viz. KsTpe = Ep e = (ft 2 /2m)(37r 2 n e o) 2 ^ 3 ), and 
as a result the ion Fermi energy should decrease to conserve the energy law. The decrease 
of the ion Fermi energy decreasing the nonlinearity of the system and hence the height of 
the soliton pulse shrinks. 

It may be interesting to note that for certain values of plasma parameters the solitary 
pulse convert to an explosive/blowup excitation as shown in Fig. 3. The blowup excitation 
indicates that an instability in the system can produce due to the effect of the nonlinearity 
(which in our case depends on the positron concentration fi p and the quantum statistical 
effects a). On the other hand, the magnitude of some quantities (e.g. temperature, pressure, 



density, etc.) leads to prejudice the balance between the dispersion and the nonlinearity [31 ]. 
Therefore, the amplitude may increase to very high values, which gives rise to increasing 
the electric potential and then accelerate the moving particles. 

It is important to notice that Eq. (23) is an explosive/blowup solution, i.e. the potential 
ip infinitely grows at a finite point (for any fixed X, Y, Z — > Xq, Yq, Zq), there exist an £o at 
which the solution (23) blowup and thereby we regard the latter as an explosive solution as 
depicted in Fig. 4. 

B. Shock/Double Layer Excitation 



For the shock/double layer solution [32j, the boundary condition <p(£) — > at £ — > oo 
must satisfy. Applying the last boundary condition into Eq. (21), we obtain the double 
layer solution as 

<p = ip m [1 + tanh (W D £)] , (24) 
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where the amplitude of the double layers is ip m = —A /2B Q , the width is Wd = 
a/— 247.E?o/Aq. Here $(= — Aq/6-Bo) is the shock wave speed. Notice that B < has 
to be fulfilled, in order for making the width Wd real. The numerical analysis in Fig. 5, 
however, shows that for small positron concentration \i v the dominant situation corresponds 
to B < 0, so the double layers may exist. For large positron concentration /i p , double layers 
cannot occur, since B > 0. Typically, we have used the plasma density value for white 



dwarf 



ll| via Tito = 2 x 10 32 cm 3 and assume that L z = 0.2, which leads to the fact that for 



negative B (i.e., formation of double layers) the positron concentration n p0 must less than 
1.43308 x 10 31 cm -3 . Also, it noted that the narrow range of \i v [corresponding to Bq < 0] 
will not change the ion gyrofrequency Q. Generally speaking, one can also note from Eq. 
(24) that the nature of the double layer depends on the sign of Aq, i.e. for A > a positive 
double layer exists (viz (p m > 0), whereas for A < we would have a negative double layer 
(ip m < 0). For white dwarf plasma parameters, it is found that A is usually greater than 
zero and then only positive double layers can exist. 

Equation (24) describes the double layer potential, which has a well-know profile (cf. Fig. 
6). This profile may change due to vary of physical parameters. The dependence of double 
layer characteristics on the positron concentration fi p [which determines Tp ep , H e and a 
through the charge-neutrality condition] is depicted in Fig. 7. It is obvious that an increase 
in the positron concentration \x v shrinks the double layers width but the amplitude increases 
by increasing \x v . 



It important to note here that in Ref. |33f ] . the soliton excitation in e-p-i 
magnetoplasma was investigated but the present work investigates soliton, shock 
and explosive excitations in e-p-i magnetoplasma. Therefore, the present model 



studies another two nonlinear structures, which did not discuss in Ref. 33]. 
Also, in Ref. 27]], the authors used the extended Conte's truncation method 



to obtain the solitary, explosive, and periodic solutions of the QZK equation. 
Note that this method gives solitary and explosive excitations described by 
equation (25) and periodic excitation described by equation (26). Thus, the 
extended Conte's truncation method cannot predict the shock formation, which 
may arise due to the presence of weakly double layers. In the present work, 
we have used generalized expansion method. The later succeeded to describe 
soliton, explosive, as well as shock excitations. Therefore, the present method 
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can be considered as a powerful tool to deal with more general nonlinear partial 
differential equations. 



V. SUMMARY 

To summarize, we have presented the properties of the nonlinear structures QIA waves 
in a very dense Fermi plasma, composed of the electrons, positrons and positive ions. By 
employing the reductive perturbation method, an extended quantum ZK equation is derived. 
The latter has been solved using the generalized expansion method to obtain a set of ana- 
lytical solutions, which reflects the possibility of propagation of various nonlinear structures 
(viz. explosive, soliton and shock pulses). We have numerically examined the effects of the 
positron concentration (which changes the quantum statistics and quantum diffraction pa- 
rameters through the charge-neutrality condition) on the electrostatic potential excitations, 
by varying relevant physical parameters. It is found that the amplitudes and widths of the 
nonlinear structures are significantly affected by the positron concentration, quantum statis- 
tical, and quantum tunneling effects. Also, for certain plasma parameters the solitary pulse 
transforms to blowup pulse. Finally, we stress that this investigation should be useful for 
understanding the features of the nonlinear structures QIA waves in an electron-positron-ion 
plasma, such as those in the superdense white dwarfs and in the intense laser-solid matter 
interaction experiments. 

Appendix: The necessary condition to neglect the annihilation 
process 

To neglect the annihilation process, the following inequality must satisfy 

— « T ann , (Al) 
where (l/u pe ) is the electron plasma period and T ann is the annihilation time. For nonrela- 



tivistic plasma, the time of annihilation reads 34 j 



3aTfi e c 



6 



1 + 60 



(A2) 



where <jt (= 6.65 x 10 cm ) is the cross section and G(= KbT /mc ) is the temperature 



range, which satisfy the inequality 



34] 
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a 2 < 6 < 1, (A3) 

where a (= 7.2974 x 10~ 3 ) is the Fine-structure constant. Equation (A3) can be rewritten 
in terms of temperature as 

3 x 10 5 < T (K) < 5.9 x 10 9 (A4) 
Inserting Eq. (A2) into (Al), we obtain 

6 > 2.66 x 10- 19 n 1/2 . (A5) 

Using Eq. (A3) and (A5), one can calculate the range of the density where the annihilation 
can be ignored 

3.9 x 10 28 < n e (cm" 3 ) < 1.4 x 10 37 . (A6) 

The quantum effects become important for certain values of density (n etP ) and temperature 
(T 6jP ). The quantum condition n e , p \\ ^ 1 specifies the temperature-density relation, where 
the quantum effects become important as 

T 6iP < 3.2 x lO -11 /^ 3 . (A7) 

Using Eq. (A6) with (A7), one can calculate the range of temperature in quantum plasma 

as 

3.6 x 10 8 < T e>p (K) < 1.8 x 10 14 . (A8) 
t is clear that the range for neglecting annihilation is well satisfied for white dwarf [see Ref. 



ll|]]. Therefore, the present model can be applicable to the dense white dwarf. 
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Figure Captions 

Figure 1 (color online): 

Three-dimensional profile of the solitary pulse [given by Eq. (22)]. A positive solitary 
pulse for k = -1, fx p = 0.8304, a = 1.693, Q = 0.05, H e = 0.03, T = 0, Y = 0.1, L x = 0.01, 
and L z = 0.1. 

Figure 2 (color online): 

Two-dimensional profile of the solitary pulse [given by Eq. (22)]. A positive solitary 
pulse for k = —1. For curve A, fi p = 0.5, a = 2.08, Q = 0.01, and H e = 0.0075, for curve 
B, fi p = 0.75, a = 1.75, Q = 0.0102, and H e = 0.0068 and for curve C, fx p = 1, a = 1.587, 
Q = 0.010291, and H e = 0.00624. Also, we have used T = 0, X = Y = 0.1, L x = 0.01, and 
L 2 = 0.1. 

Figure 3 (color online): 

Three-dimensional profile of the explosive/blowup pulse [given by Eq. (22)]. A positive 
explosive pulse for k = -1, fi p = 0.6, a = 1.9, Q = 0.0102, H e = 0.0072, T = 0, Y = 0.1, 
L x = 0.01, and L z = 0.1. 

Figure 4: 

Three-dimensional profile of the explosive/blowup pulse [given by Eq. (23)], for ji p = 
0.0525, a = 7.37, = 0.0257, H e = 0.0081, T = 0, Y = 0.1, L x = 0.01, and L z = 0.2. 
Figure 5: 

The nonlinear coefficient B is depicted against the positron density n p0 for n i0 = 2 x 10 32 
cm -3 and Lz = 0.2. Recall that for n p0 < 1.43308 x 10 31 cm -3 the nonlinear coefficient B < 
and then a shock pulse can propagate. 

Figure 6: 

Three-dimensional profile of the shock pulse [given by Eq. (24)], for fi p = 0.05, a = 7.6, 
tt = 0.03, H e = 0.008, T = 0, Y = 0.1, L x = 0.01, and L z = 0.2. 
Figure 7: 

Two-dimensional profile of the shock pulse [given by Eq. (24)]. For curve A, fi p = 0.05, 
a = 7.6, VI = 0.03, H e = 0.008 and for curve B, fi p = 0.052, a = 7.37, = 0.03, and 
H e = 0.0081. Here, T = 0, X = Y = 0.1, L x = 0.1, and L z = 0.2. Recall that the narrow 
range of ji p will not affect on the ion-gyro frequency Q. 
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